Spectral theory and distributional Borel summability 
for the quantum Henon-Heiles model 



Emanuela Caliceti 

Dipartimento di Matematica, Universita di Bologna 40127 Bologna, Italy 

and INFN, sezione Bologna 
(caliceti@dm.unibo.it) 



Abstract 

The Borel summability in the distributional sense is established of the divergent 
perturbation theory for the ground state resonance of the quantum Henon-Heiles 
model. 

1 Introduction and statement of the results 

A standard model for transition to chaos is given by the Henon-Heiles Hamiltonian 
(see e.g. [16], [12], [18]) denned by 

H(P)=pl+pl + q l + £ + p(q 2 1 q 2 -~qi S j, G R. (1.1) 

The quantum counterpart of (1.1) (see e.g. [15], [2], [1], [11], [3], [14]) is represented 
by the Schrodinger operator in L 2 (R 2 ) formally given by 

H(J3) = -A + x 2 + (3 (x\x 2 - := H(0) + @V, (1.2) 

where x = (xi,x 2 ) G R 2 , x 2 = \x\ 2 = x\ + x 2 , A is the 2— dimensional Laplace oper- 
ator, V is the multiplication operator by the function V(x) I and H(0) 
is the operator corresponding to the harmonic oscillator. Its spectral properties have 
been extensively investigated numerically (see e.g. [13]). However the mathematical 
analysis of the problem is somewhat tricky. The purpose of this paper is to solve one 
of the mathematical problems involved, namely the meaning of perturbation theory, as 
announced in [5], in order to obtain results analogous to those obtained in [7] and [4] 



for the one dimensional odd anharmonic oscillator. Indeed, as in the one dimensional 
case, the minimal operator generated by (1.2), with C£°(R 2 ) as domain, is not essen- 
tially selfadjoint; in fact, it has infinitely many selfadjoint extensions, none of which 
with discrete spectrum. Therefore the numerically observed eigenvalues have to be 
interpreted as (real part of) resonances [19], which are complex objects: the real part 
is the location, and the imaginary part the width. Another singularity in this problem 
is due to the fact that the Rayleigh-Schrodinger perturbation expansion (from now on 
denoted RSPE) near the lowest unperturbed eigenvalue Eq = 2 not only diverges but 
has real coefficients with constant signs. As is well known, this prevents the series from 
being Borel summable in the ordinary sense. The aim of this paper is then to prove 
that, as in the Stark effect (see [9]), the perturbation series near E is Borel summable 
in the distributional sense to the real part of the resonance. On the other hand, if f3 is 
taken to be purely imaginary, then H(/3) shows a less singular bahavior (see [19] and 
[6]). In fact H{(3) is closable on Cq°(R 2 ) and, if D(H((3)) denotes the domain of its 
closure, then H(0) is PT-symmetric, i.e. 

(H(/3)u)(-x) = H(py^Px), Vu G D{H{(3)). 

Furthermore the RSPE near any eigenvalue has real coefficients with alternating signs 
and is Borel summable to the corresponding (real) eigenvalue of H(f3) (see [19]). Thus, 
the natural way of dealing with H((3) is to start with j3 complex and then look for a 
continuation to (5 real. More precisely we will prove the following. 

Theorem 1.1 (a) H((3) represents an analytic family of type A of closed operators 
with compact resolvents, with domain 

D{H{(3)) = D(H(0)) n D(V), for < ar g/ 3 < tt. 

(b) Let E = 2 denote the ground state energy level of H(0). Then for any 5 > 
there is B(8) > such that for \/3\ < B(5), < arg/3 < ir, H((3) has exactly one 
eigenvalue E(f3) near E , which admits an analytic continuation across the real 
axis to the sector 

S 5 := j/3 G C : < \(3\ < B(8), ~ + 5 < arg/3 < - 5 J . (1.3) 

Moreover, lim^o E(/3) = E 

In [19] weaker results concerning the spectral properties of H(f3) are obtained, which 
are not sufficient to establish our main result, stated in Theorem 1.2 below, i.e. the 
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distributional Borel summability (from now on denoted DBS) of the RSPE around Eq. 
More precisely in order to specify the domain of H(/3) and guarantee the compacteness 
of its resolvents (as stated in Theorem l.l-(a)), we need a quadratic estimate on H(/3), 
proved in Appendix B. We can now state the main result of this paper. 

Theorem 1.2 Let (5 e R. Then 

(a) the RSPE near Eq is Borel-Leroy summable of order | in the ordinary sense to 
E{($) for < arg/3 < -n and in the distributional sense to $tE(ft) for (3 e R, \(3\ 
suitably small; 

(b) UE{(3) = &E(-P), $fE(/3) = -9f E(-/3), for(3e R. 

The notion of distributional Borel-Leroy summability and the corresponding criterion 
are recalled in Appendix A for the convenience of the reader. 

Remark 1.3 In order to extend the result stated in Theorem 1.2 to any unperturbed 
eigenvalue E t — 2(1 + 1), I — 1, 2, ... of H(0), one needs to extend the definition (and 
corresponding criterion) of distributional Borel summability to the degenerate 
Hunziker and Pillet did in [17] for the notion of ordinary Borel summability. We plan 
to do this in a forthcoming paper. 

The proof of Theorem 1.1 is obtained in Section 2, following [4] where analogous re- 
sults were obtained for the one dimensional add anharmonic oscillator. We will include 
most of the details in order to make the paper self-contained. The proof of Theorem 
1.2 requires the verification of the analogue of the Nevanlinna criterion for the DBS, as 
stated in Appendix A and proved in [8]. In particular we need to extend the analyticity 
of E(f3) to a suitable Nevanlinna disc. This result is achieved in Section 3 by using a re- 
finement of the Hunziker- Vock stability technique introduced in [9], [10] and [4] where 
analogous results are obtained for the resonances of the Stark effect, of double- well 
oscillators and of the 1— dimensional odd anharmonic oscillator, respectively. However 
such technique cannot be directly applied to the present 2— dimensional problem. In- 
deed, the method is based on a control of the numerical range of H(/3), whose distance 
from any complex number z o~(H(0)) must be bounded from below by a positive 
constant as \x\ — > oo and f3 — > 0. It is possible, however, to overcome this difficulty 
by passing to polar coordinates, where the angular coordinate can be easily dealt with 
and the problem essentially reduces to a 1— dimensional one. 
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2 The operator H(/3) for ^ (3 ^ and the analytic 
continuation of the eigenvalues 

In order to prove Theorem 1.1 we need some preliminary results. 
Lemma 2.1 Let f3 = \f3\e ta with a G]0, 7r[ and Q be a compact subset of 

{7 G C \ {0} : |7| 2 > 4|/3| sin a, -tt + a < arg (3 < a}. 
Then there exist a, b > such that 

\\Au\\ 2 + ItHI^H 2 + |/5| 2 ||^«|| 2 < a||(-A + ix 2 + (3V)u\\ 2 + 6||«|| 2 (2.1) 

\/u G C^°(R 2 ), 7 G fi, < \(3\ < 1, a and b independent of 7 m Q and a in a closed 
interval contained in ]0, 7r[. 

Proof. See Appendix B. 

Corollary 2.2 Let 7, /3 G C satisfy the conditions of Lemma 2.1. Then the operator 
T( 7 ,/3) denned in L 2 (R 2 ) fry 

T( 7 , = -Am + 1X 2 u + /3Vu, Vm G D(T(rf, (3)) 

on the domain D(T(j,{3)) = D(H(0)) fl D(V) is closed and has C£°(R 2 ) as a core. 

Lemma 2.3 Let 7, /3 G C satisfy the conditions of Lemma 2.1, i.e. a = arg/? G 
]0, 7r[, — 7r + ct < arg7 < a, \^\ 2 > A\/3\ sin a. 
Then there exists £ > snc/i £/ia£ 

(Sft [e-^^ 2 ' < «, T( 7 , /3)m >]><«, - Am > G D(T(>y, /?)). (2.2) 

Proof. It is enough to prove (2.2) for -u G C^°(R 2 ). We have 

SR [e"*^/ 2 ) < u , (-A + 7 x 2 + 0V)u >] 

= cos (a — —^j < u, —A-u > +I7I cos ^— — a + arg 7^ < u, x 2 u > 

= sin a < u, —Aw > +\*y \ sin(o; — arg 7) < u,x 2 u >> sin a < u, —Au > 

since sin(o; — arg 7) > because < a — arg 7 < n by assumption. Now, since 
< a < n, the lemma is proved with £ = (sina) -1 . 

Corollary 2.4 TTie numerical range ofT{j,f3) is contained in the half-plane {zgC: 
— 7r + a < arg z < a} . 



Set H{(3) := T(l, f3), for a = arg (3 G]0, ir[. By the above results we have the following 
theorem, which corresponds to Theorem 1.1- (a). 

Theorem 2.5 H(f3) represents an analytic family of type A of operators with compact 
resolvents, with D(H(f3)) = D(H(0)) n D(V), for < arg/3 < vr. 

We will prove that the (discrete) spectrum of H(/3) is non-empty (see also [19] ). In 
order to prove Theorem 1.2, i.e. the DBS of the RSPE near Eq we need to prove 
that it is stable with respect to the family H(j3), as (3 — > 0,3/3 > and that the 
corresponding eigenvalue E((3) of H(/3) can be analytically continued to a wider sector 
than < arg (3 < ir. To this end we start by making use of standard dilation analyticity 
techniques, i.e. we introduce the operator 

H((3, 9) := -e- 2e A + e 2 V + (3e se V(x) := e- 2e K((3, 9) (2.3) 

which, for 9 G R, is unitarily equivalent to H((3), ^s(3 > 0, via the dilation operator 

(U(9)u)(x) = e 9 u(e 9 x), G L 2 (R 2 ) 

First of all notice that K(f3, 9) = -A + e Ad x 2 + f3e 5e V(x) corresponds to T(j, (3') with 
7 = e 4e ,/3' = f3e 56 '. From now on we will assume \$t9\ < 1 and \(3\ < ^ so that the 
conditions \(3'\ < 1 and \^\ 2 > 4|/3'| sin(arg/3') required in Lemma 2.1 are automatically 
satisfied. 

Next observe that the further conditions on (3' and 7 so far required: 

-vr + arg (3' < arg 7 < arg (3' ^ 4 ^ 

< arg (3' < 7r ^ ' ' 

are equivalent to the following conditions on (3 and 9: 

O<arg/3 + 30<7r , , 

< arg/3 + 5^ < vr 1 J 

In complete analogy with Theorem 2.5 of [4] we can now prove the following 

Theorem 2.6 Let s = arg (3,t = ^s9. Then H((3, 9) is a holomorphic family of type A 
of closed operators on D(H((3, 9)) = D(H ) n D(V) with compact resolvents for (3 and 
9 such that s and t vary in the parallelogram P of the (s,t)— plane defined by 

P = {(s, t) G R 2 : < t + s < 7T, < 5t + s < vr} . 
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Remark 2.7 1. By Corollary 2.4, the numerical range of K(/3,8) is contained in 
the half-plane — n + a < argz < a with a = arg/3 + 530; thus, if 0) has 
numerical range contained in the half-plane 

II = {z G C : -7T + arg/3 + 330 < arg^ < ar g/ 3 + 330} . (2.6) 

Moreover the (discrete) spectrum of H(j3, 9) is contained in II and ^ II, || (z — 
H(f3,9))- 1 \\ < dist(z,II)" 1 . Finally, the analyticity of H(f3,9) in the region 
defined by P allows f3 to be extended to the sector S — {(5 : < \/3\ < j3 , — | < 
arg/3 < fvr} 

2. If we start from the operator if(/3) with 3/3 < 0, analogous results can be ob- 
tained for the operator family if (/?, 9) for (3 and 9 satisfying the conditions: 

-7T < 30 + arg/3 < , , 

-7T < 530 + arg/3 < 1 ' 

Furthermore H{(3, 9)* = H{fi, 9). 

By standard dilation analyticity arguments the eigenvalues E\ — 2(1 + 1), I — 0, 1, . . ., 
of if (0,0) = -e- 2e A + e 26 x 2 ,D(H(0,9)) = D(H(0)), are independet of 9 for -f < 
30 < ^. By an argument similar to that used to prove Theorem 2.7 in [4] (see also 
[19]) we can prove the stability in the sense of Kato of each eigenvalue E\ with respect 
to the family {H({3, 9) : \/3\ > 0}, (5 and 9 in the region defined by P. More precisely 
we can state the following 

Theorem 2.8 Let (3 and 9 satisfy conditions (2.5). We have 

(a) if\£ a(H(0,9)), then XeV, where 

V = {z G C : z £ a(H(P, 9)) and(z-H((3, 0))" 1 is uniformly bounded as \(3\ -> 0} 

(b) if A G cr(if (0, 9)) = {2(1 + 1) : 1 = 0,1,...}, then A is stable with respect to the 
family H(j3,9), i.e. if r > is sufficiently small, so that the only eigenvalue of 
if (0,0) enclosed in F r — {z G C : \z — A| = r} is \, then there is B > such 
that for \(3\ < B, dimP(/3, 0) = dimP(O,0) ; w/iere 

0) = (2m)- 1 £ (z - if 9)Y 1 dz 

is the spectral projection of H((3, 9) corresponding to the part of the spectrum 
enclosed in Y r C C \ a(H(0, 9)). 
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Remark 2.9 It can be immediately checked that all the results so far obtained, in 
particular the analyticity of the family H(/3, 9) and the stability of the eigenvalues of 
the harmonic oscillator H(0, 9) with respect to H((3, 9) as \/3\ — > 0, hold unifomly in (3 
and 9 such that (arg/3, 9) varies in any compact subset of P. 

Now we specialize the result obtained in Theorem 2.8 to the ground state energy level 
E = 2 of H(0). More precisely, for any 5 > there exists B(5) > such that for 
\(3\ < B(5), — \ + 5 < arg/5 < | n — 8, H(/3, 9) has one and only one eigenvalue E((3), 
independent of 9 if (arg/3, Q6) G P, which converges to E as — > 0. 
By Theorem 2.6, is analytic in the sector 

S 6 = \j3 G C : < < + 5 < arg/5 < ^tt - <j| (2.8) 

and is an eigenvalue of H((3) for < 3/3 < n. For future reference we state this result 
in the following 

Theorem 2.10 For any 5 > 0, there is B(5) > such that for < B(5), < 
arg/5 < 7T 7 i?(/3) has exactly one eigenvalue E(/3) near E , which admits an analytic 

continuation across the real axis to the sector S$- Moreover lim^o E(f3) = E Q . 

(3es s 

3 Analyticity of E(/3) in a Nevanlinna disc and DBS 

The basic analyticity result needed to establish the DBS of the RSPE near E for (3 > 
is obtained in the following theorem. 

Theorem 3.1 There exists R > such that the eigenvalue E((3) of H(f3) = —A + 
x 2 + (3{x\x 2 — \x\) near E for \/3\ small is analytic in the Nevanlinna disc Cr = {(3 : 
> R- 1 } in the (3 2 -plane. 

Remark 3.2 (I) The sector S{8) given by (2.8) can be rewritten in terms of the 
variable j3 2 as: 

S{5) = j/3 : < \/3\ < B{5), ~ + 25 < arg/5 2 < | + 2tt - 25 j . (3.1) 

(II) The function E{(3), analytic in any sector S(5) and for which we want to prove 
analyticity in Cr, represents an eigenvalue of H((3,9) if (arg/3, Q6) G P. In 
particular, for — | < arg/3 < we can choose the path inside P given by the 
straight line of equation 

%9 = - -arg/3 + -. 
3 6 
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Then, if we set 



argp = 1 — , i.e. argp = h e, e — ► IT 

sp 4 2 2 



(3.2) 



we obtain S6> = | — |, and the operator if 6 1 ) takes the form 
A{p) = -e 2 e-^-^A + OfeW-Vx 2 + ip6lV{x) 
with p = \(3\ and O = e 3 ^. 
(Ill) For /9 = pe iarg/3 and arg/3 = — | + |, the boundary of (7r has the equation 



(3.3) 



sine = — . 
R 



(3.4) 



(IV) Since the disc Cr is the union of the boundaries of discs of smaller radius, the 
proof of Theorem 3.1 reduces to a stability argument for the eigenvalue E with 
respect to the family {A p } p>0 as p — > + . In view of Remark 3.2 (IV) it is 
convenient to move to polar coordinates (r, ip) as follows: 



I x\ — r cos if 
| x 2 = r sin if 

Then A(p) is equivalent to the operator Ai(p) formally given by 



(3.5) 



d 



1 d 2 



1 



)_ 02gi(f-f ) r 2 



+ ipO^r 3 ^cos 2 sin if — ^ sin 3 ip 



(3.6) 



in L 2 (R + x T 1 ) with suitable conditions at the origin. 

More precisely A 1 (p) is the operator generated by the quadratic form 



a p [u\ 



oo r2ir 

o Jo 



M 2 + ^MI 2 -^M 2 



+ 



i 



#oe i( 2 a) r 2 + ^p6)3 r 3 (cos 2 v9sinv9 - - sin 3 <p 



|u| 2 jdrdfl 



(3.7) 



defined on the maximal domain with the condition at the origin u(r, if) ~ r 1 / 2 
uniformly in ip. Here we adopt the notation p r = —i-j^^p^ = —i-§^- Now, as in 
[9], let U be a transformation in the set of L 2 functions which are translation 
analytic in some suitable strip |3r| < 5q, defined by: 



{m){r,d)=i\rf/ 2 ^{r)^) 



(3.8) 
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where, setting r = y for a suitable do > 0, we define £(r) G Cq°(R + ) so that, 
for ?7o e]0, 1[ fixed, 

^(r) =r-2^ [l-(l + r 3 )- 1 / 6 ], < r < r (3.9) 

£(r) = r, r>r + 77o (3.10) 

Setting /(r) = |7^y and = e~ % ^~i\ the transformed operator UA(p)U~ 1 is 
given by 



Prf 2 Pr + \(f 2 )"-^pl-^ 



+ uj 1 6l£ 1 + ip6^ 3 (cos 2 y? sin y? - ^ sin 3 y?) . 
The quadratic form which ganerates H p is defined by: 



(3.11) 



h p [u\ 



roc r'27v 






u6f 


lo Jo 



f\PrU\ 2 + -^\P^U 



+ 



4 

33 f 3/ „ 2 



\u\ 2 \drd<p (3.12) 



on the maximal domain with the condition at the origin u(r, ip) ~ r 1 / 2 , uniformly 
in (p. The limit in the strong resolvent sense as p — > + of H p is defined by 



1 



(3.13) 



where / = (Co) -1 and £ (0 is defined by (3.9) Vr > 0. From Remark 3.2-(IV) 
it follows that in order to prove Theorem 3.1 it suffices to prove the stability of 
the eigenvalue E of H with respect to the family H p as p — > + . As in [9] this 
result is achieved by means of some preliminary results. 



Lemma 3.3 Let sine 
ai, a2, c, Ci, C2 such that 

/•2n rrg 



R 



There exists p > 0, n e N and positive real constants 



m p [u] > ai JJ £ | sin (0 [1 - r 4 (l + r 3 )- 7/3 ] + 7? r 2 (l + r 3 )- 7/6 ||p rM | 2 rfr^ 



+a 2 



2-7T /"OO 



Jr 



sin ( - ) \p r u\ 2 drd(p — c\\u\ 



(3.14) 
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Vp e]0,po] 5 Vit G D(h p ). A similar estimate holds with h p replaced by ho (r = +00, e = 
0). 

Moreover, Vu e D(h p ) such that supp w C (n, +00) 



9*M«] ^ (ciiT 1 -c 2 )|H| 2 
Vn > n and c±, c 2 independent of R, Vp G [0,p ]. 



(3.15) 



Proof. We have 

f27T POO \ 



f\Pru\ 2 - —\u\ 2 \drdip 
Jo I J 



dip f 
Jo Jo 



bi2 1 1 1 o 



(3.16) 



Since the function ^ — |z is bounded we obtain 



3? 



I Jo " y 1 ^ 2 - 4e |M|2 ^| 



/•27T /'(DO 

> / / ®("/ 2 ) 
JO JO 



r2w pro 



drdip — (const.) || m I 

2w roo 



rzn H rZir pro rzir roo 

= / / (■■■)drdip+ / / (■■■)drdip+ / / (• • -)drd<p - (const.) || 

Jo JO JO Jl JO Jr 



'II- 



(3.17) 



Let us now denote Ii, I 2 and J 3 the first, second and third integral respectively, in the 
right hand side of (3.17). For r < r we have 



> ^{l-r^l + r 3 )- 7 / 3 } 
3/ 2 > 2-V 2 (l + ^ 3 )" 7/6 - 



(3.18) 



Hence, 



^(co/ 2 ) > A' 1 sin (|) [1 - r 4 (l + r 3 )" 7 / 3 ] + ^ cos (|) r 2 (l + r 3 )" 7 / 6 . (3.19) 



Now, by Sobolev's inequality I\ > 0; moreover 

f 2-7T /"OO 



h > J J sin (^-J |p r -u| 2 <ir<i<p — (const.) ||w|| 2 , 
I 2 > a, fj £ {sin Q) [1 - r 4 (l + r 3 )- 
+r]o^ 2 (l + r 3 )~ 7/6 j \p r u\ 2 drdip . 



.3 \ -7/3] 



(3.20) 



(3.21) 
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Let us now estimate the remaining terms appearing in $th p [u]. We have 

2-7T pOO _ /•27T pOO J*2 4^^(r) 



9? 

where 



L L Z-Wri^l I J ^± M ^ (3.22) 



77o[l - (1 + r 3 ) 1/6 ], for r < r 



(r)= w ■ j i, — - --u (3>23) 
/v y [ for r > r + r? • 

By choosing r] suitably small (e.g. r] < 2~ 6 / 5 ) we have (r 2 — Arf){r 2 + r? 2 ) 2 > and 
therefore (3.22) is non negative. Next notice that (f 2 )" is bounded; thus, 

& / ^^(f)"\u\ 2 drdp> -(const.) |M| 2 . (3.24) 
Jo Jo 4 

We can now estimate the potential term V p (r, ip) := a~ 1 9 2 ,£ 2 + i p9^ 3 (cos 2 p sin p — 
| sin 3 ip). We have 

K(V p (r,p)) = sin (0 (r 2 - V (r) 2 )9 2 + 4cos (|) r^r)^ 2 

+6pr 2 #Q?7(r) ^cos 2 y9 sin (/? — ^sin 3 ^ . (3.25) 
Fixing R > 0, for any 6 > there exists k > such that for r e (0, 6): 

MV p (r,p)>-k. (3.26) 

Moreover: 



(I) For r > r + ??o? recalling that r = a /p, we have 

R r °~~R 



^V p (r,p)>k4r 2 Q = ^ (3.27) 



for a suitable constant k x > 0. 

(II) Finally, for 6 < r < r + r] (6 > 0, independent of p > 0) we have < ??(6) < 
T](r) < r] . Thus, 

5 * "2„/„\ o„2_ /i3 



3^(r,<p) > -A; 2 + 4cos(^-Jr^(r)-8r 2 p^ 3 

> -fc 2 + 2rr](b)6 2 - 8r] r 2 p6 3 (3.28) 

for some k 2 > 0, if p > is sufficiently small so that cos(|) > \- By suitably 
choosing a > (e.g. a = r](b)/8r]o9 , 6=1) the term 2ri](b) — 8pi] r 2 attains its 
maximum at r = r = M and its minimum at r = b in (6, r ). Thus, for r G (6, r ) 
we have 

$V p (r, p) > -k 2 + 2b 2 r](b)6 2 - 8pr] b 2 6 3 > -k 2 , (3.29) 
and this concludes the proof of (3.14). 
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As for (3.15) notice that the kinetic part is > — (const.) \\u\\ 2 , and for the potential 
part we have: 

(F) For r > r + r] we can repeat the argument used in (I). 

(IF) For n < n < r < r we proceed as in (II) with n > b {b = l),i](r) > rj(b); thus, 

®V p (r, <p) > -k 2 + 2r V (b)6 2 - 8r] r 2 P e 3 . (3.30) 
Again 2rrj{b)6 2 l — 8^oP r2 ^o attains its maximum in (n, r ) at r and its minimum 



at r = n. Hence 



®V p (r,(p) > -k 2 + 2nri(b)9 2 - 8n 2 pri 9 3 

k 5 

> k 3 n — k A > — — k%. (3.31) 
rt 

This concludes the proof of the Lemma. 

Corollary 3.4 Let x„(r) = x( r / n ), X £ C°°(R + ), x( r ) — 1 f or r < 1> x( r ) — if 
r > |. T/ien £/iere exists C3 > snc/i £/icrf 

||[#„XnHI < c 3 n- 1/4 (||# p u|| + ||«||) Vm e < p < Po . (3.32) 

Proof. It is enough to prove (3.32) for p > 0, since for p = the argument is simpler. 
For simplicity we set ||w|| = 1. Let 72n(?") be the characteristic function of the interval 
[1, 2n] in the r variable. Then we have: 

[H P , Xn] = -e, 2 ^ [2in- l f 2 X > (£) + n- 2 X " (£) + 2//Y (£) } ■ (3.33) 



Notice that the term — ^|§-p^ in H p gives no contribution to the commutator [H p , x n ] 
since x( r ) does not depend on <p. Now: 



{/ /-27T r2n \ V 2 



< c 5 n 1 



2n 



'l + 8// :1 ) 7 " 2 , /■*'< , , „ „. . .„_ 7 , R l 1/2 + 1 ' 



y dip J dr\p r u\ 2 r] r 2 {l + r 3 )- 7/6 





1/2 1 


\ 2 dr 


+ 1 l 



< c 7 n~ 1/4 {[K < H p u,u> +c 8 ] 1/2 + 1} (3.34) 

where the last inequality follows from Lemma 3.3. Now, taking c$ = eg + 1 with 
3? < H p u, u > +c 9 > we obtain: 

II [H p , Xn\u\\ < Cion- 1 / 4 !^ < H p u, u > +c n } (3.35) 

where c\\ = c$ + 1. This concludes the proof. 
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Remark 3.5 A similar argument can be used to obtain the analogous estimate for the 
adjoint operator H*\ 

\\[H* p ,Xn}u\\ < cgn-^dl^ull + HI) Vu G D(H*), < p < p . (3.36) 

Proposition 3.6 Let M n = 1 — where Xn is defined as in Corollary 3.4- 
U 

d n (\,p)=m£{\\(\-H p )M n u\\ : ||M n u|| = 1, u G D(H P )} 
then VA G C, 3R, no, po, S > such that 

d n (X, p) > 5 > 0, Vn > n , Wp < p . (3.37) 

Proof. First of all notice that 

d n (X,p)>dist(X,E n (p)) (3.38) 

where 

E n (p) = {< M n u, H p M n u >: u G D(H p ), \\M n u\\ = 1} 
Moreover, by (3.15) we have: 

3? < M n p, H p M n u >> CxR' 1 - c 2 . (3.39) 

Now, (3.37) follows from (3.39): since C\ and c 2 are independent of R we can take 
R > suitably small so that (3.37) is satisfied with S — ^ — c 2 — |A|. 

Lemma 3.7 Lei p m > andu m G D(H Pm ) be two sequences such that p m + , ||if Pm u 
is bounded and \\u m \\ = 1, u m A 0. T/ien £/ie sequences p m i n ), M m u m ( n ) satisfy the same 
properties for suitable m = m{n), if R> is chosen sufficiently small. 

Proof. By Corollary 3.4 the boundedness of ||if Pm « m || implies that of \\H Pm{n) M n u m ^ n )\\. 
Thus it is enough to prove that 

lim ||Xn«m(n)|| = 0, Vn . (3.40) 
To prove (3.40), let H' p = ^ l H p and A G C \ a(H' Q ) be fixed. Then 

\\XnU m \\ 2 < c(\\XnK(H' - #;>m|| 2 + \\XnR'o(H' Pm ~ X ) U ^) " ( 3 ' 41 ) 

The second term in the right hand side of (3.41) tends to zero as n — > oo, because 
i? := (ifg — A) -1 is compact and {H' — \)u m A 0. The first term in the right hand 
side of (3.41) can be bounded, up to a constant factor, by 

\\R' Xn{H' - H' pm )u rn \\ 2 + \\[R' , XnKK - H' pm )u m \\ 2 . (3.42) 



13 



Now, the first term in (3.42) can be bounded as follows 
\\R' oX n(H' - H' p Ju m \\ 2 < c\\R' \\ 2 dp \(H' -H' p Ju m \ 2 dr 

JO J 

< c\\R> \\ 2 dtp (\^r)\ 2 \l-e-^\+p m \i{r)\^f\u m \ 2 dr 

Jo Jo 

< c'(e m n 2 + p 2 m n 3 ) 2 \\u m \\ 2 (3.43) 

where sin(£ TO ) = p^/R by hypothesis. In the second inequality we have used the fact 
that, for r G (0, 2n), £(r) = £o( r ) = r — 2ir)(r), r)(r) = i] [l — (1 + r 3 ) -1 / 6 ] < r] . Since 
(3.43) tends to zero as n — > oo, let us estimate the second term in (3.42): 

||[i^,X»](^o-^)«mll < KWx»TO^o-^)«mll <^n _1/4 - (3-44) 

Indeed, the operator R' [Hq, Xn] = ([Xn, (Hq)*](R' )*)* , when applied to the bounded 
sequence R' {H' — H' )u m satisfies the inequality (3.44) by (3.36). 

Proof of Theorem 3.1. Since lim p _^ + H p u = H u, lim p ^ + H*u = HqU, Vm G C£°(R + x 
T 1 ), we can use Corollary 3.4, Proposition 3.6 and Lemma 3.7 in order to apply The- 
orem A.l of [10] which provides the following stability result: 

(1) if A ^ cr(H ) then (A — Hp)" 1 is uniformly bounded as p — > + ; 
(ii) if A G cr(H ) then A is stable with respect to the family {H p } p>0 . 

Now the proof of the theorem is a consequence of Remark 3.2-(IV). 

Proof of Theorem 1.2. Taking into account the result obtained in Theorem 3.1, the 
proofs of Theorems 3.13 and 1.3 of [4] can now be taken over directly without change 
in order to prove (a) and (b) respectively. 

Remark 3.8 (1) For (3 G C R ,E(f3) and E{J3) are the so called "upper sum" and 
"lower sum" respectively of the RSPE (see Remark A. 3 below), while the distri- 
butional Borel sum is given by f((3) = \{E{(3) + E(p)) and d{(3) = E{(3) - E(0) 
is the discontinuity with zero asymptotic expansion. The result obtained in The- 
orem 1.2 can be interpreted in terms of resonances of the problem as explained 
in Remark 3.14 of [4]. 

(2) Similar results can be obtained if we now start from < 0, instead of > 0. 
We can establish a relationship between the resonance Ei(/3) obtained in this 
case and E(f3) following again [4]. Indead we have Ei(/3) = E(/3) for (5 G R. 
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A Appendix 



To make the paper self-contained, in this appendix we first recall the notion of distri- 
butional Borel-Leroy summability of order q as introduced in [8]. 

Definition A.l Let q be a rational number, (a s ) s€ N Q> sequence of real numbers and 
R > 0. We say that the formal series Y^=o a sft s is Borel- Leroy summable of order q in 
the distributional sense to f{(5) for < (3 < R if the following conditions are satisfied. 

(a) Set 

Then B(t) is holomorphic in some circle \t\ < A; moreover B(t) admits a holo- 
morphic continuation to the intersection of some neighborhood o/R + := {t G R : 
t > 0} with C+ := {t G C : > 0}. 

(b) The boundary value distribution B(t + iO) exists Vt G R + ; and the following 
representation holds: 

1 roc 

f(P) = ^J o PPm))e~ W)1/q l^j dt (A.2) 

for (3 belonging to the Nevanlinna disc of the f} 1 ^— plane Cr := {(3 : 9?/3 -1 / 9 > 
R- 1 }, where PP(B(t)) = \{B(t + iO) + B(t + iO)). 

If q = 1 the series is called Borel summable in the distributional sense to f{(3). 

Let us now recall the criterion for the distributional Borel-Leroy summability (see [8]). 
As for the ordinary Borel sum, it shows that the representation (A.2) is unique among 
all real functions admitting the prescribed formal power series expansion and fulfilling 
suitable analyticity requirements and remainder estimates. For the sake of simplicity 
we limit ourselves to the case q — 1. 

Theorem A.2 Let f(f3) be bounded and analytic in the Nevanlinna disc Cr = {(3 : 
1 > R' 1 } and let f((5) = ($(/9) - $(/3))/2, with analytic in C R and such that 

N-l 

$(/?)- a sP S < Aa(e) N N\\p\ N , ViV=l,2,... (A.3) 

s=0 

uniformly in Cr^ £ = {(5 G Cr : arg/3 > — rr/2 + e}, We > 0. Then the series 
J2"^Lo{ a s/ s\)u s is convergent for small \u\ and it admits an analytic continuation B(u) = 
B 1 (u) + B 2 (u), where Bi{u) is analytic in C\ = {u : dist(u, R + ) < d^ 1 } and B 2 {u) 
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is analytic in Cj={u: (^su > 0, > — d x ) or\u\ < d x } for some d > 0. B{u) 
satisfies 

\B(t + M)\<A , (r ]0 )- 1 e t l R (A.4) 

uniformly for t > 0, for any r) such that < i] < d^ 1 . Moreover, setting PP(B(t)) = 
(B(t + iO) + B(t + i0))/2, admits the integral representation 

f{(5) = /T 1 / PP(B(t))e-Wdt, PeC R (A.5) 
Jo 

i.e. f{(5) is the distributional Borel sum of YlT=o a sP s f or < (3 < R in the sense of 
Definition A.l. Conversely, if B{u) = J2^Lo( a s/ s\)u s is convergent for \u\ < d' 1 and 
admits the decomposition B{u) = Bi(u) + B 2 {u) with the above quoted properties, then 
the function defined by (A.5) is real-analytic in Cr and = / °° B{t-\-io)e' t ^dt 
is analytic and satisfies (A. 3) in Cr. 

Remark A.3 The function <j>(p) = (3~ l / B(t + i0)e~ t//3 dt is called "the upper sum 

Jo 

— - roo 

and 4>((3) = (3" 1 \ Bit + itye^^dt "the lower sum" of the series. It follows that, for 
Jo 

(3 > 0, f{f3) = Re0(/3). On the other hand with this method we can single out a unique 
function with zero asymptotic power series expansion, that is the "discontinuity" 

/•OO 

d((3) = /T 1 / (B(t + iO) - B{t + i0))e- t/(i dt = <j>(/3) - <j>(/3) . 
Jo 

Thus, d((3) = 2iIm0(/3), for (3 > 0. 

B Appendix 

Proof of Lemma 2.1. We shall proof the following estimate, equivalent to (2.1): 

||A«|| 2 + \a\ 2 \\x 2 u\\ + |,3| 2 ||V^|| 2 < a||(- e - ia A + ax 2 + \(3\V)u\\ 2 + b\\u\\ 2 (B.l) 

Vw G D(Hq) fl D(V), with a = ^ye~ ia varying in a compact subset of {a G C \ {0} : 
\a\ 2 > A\(3\ since, — n < arg7r < 0}. From now on we shall use the notation —A = 
pl+pl, where pj = — J^r, j — 1,2. As quadratic forms on C£°(R 2 ) (g> C^°(R 2 ) we have 

{-e ia A + ax 2 + \(3\V{x)){-e- ia A + ax 2 + \/3\V(x)) 
= {-e ia A + \(3\V{x))(-e- ia A + \/3\V(x)) + |a| 2 |a;| 4 

+$to-[(-e ia A + \I3\V{x))x 2 + x 2 (-e- ia A + \P\V(x))\ 

+i^a[(-e ia A + \P\V(x))x 2 - x 2 (-e~ ia A + \P\V(x))\ 
= — (-e- ia A + \(3\V(x) ± \a\x 2 )(-e- ia A + \(3\V(x) ± \a\x 2 ) 
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+ 1 



Ma 



[(-e ia A + \/3\V(x))(-e- ta A + \(3\V(x)) + |a| 2 |x' 41 



+i%a(-e ia Ax 2 + e~ ia x 2 A) 
Ma 



> 1- 



2|„|4i 



[(-e m A + |/3|\/(a;))(-e- m A + \(3\V(x)) + \a\ 2 \x 



+i^sa cos ct[— A, x 2 ] — ^a sin a(— Ax 2 — x 2 A) 
Ma \ , 



1 - 



+ 2^0- cos a(p 1 x 1 + XxV\ + P2X2 + X2P2) 



— Scr sin a(p 1 x 1 + x l p l + p 2 x 2 + x 2 p 2 + 2x x p 2 + 2x 2 p x ) 
(since Qa < 0, sin a > and x\p\ + x\p\ > 0) 

^ [•••] + 2^sa cos a(p 1 x 1 + XiPi + p 2 x 2 + x 2 p 2 ) 



> 1- 



Scr sin a(]9 2 x 2 + x 2 p 2 + p 2 ,:*; 2 , + ^2^2) 

[• • •] ± 2Qa\ cosa\(piXi + x x p x + P2X2 + x 2 p 2 ) 



1 - 



a 



—^sa sin a(— 4 + 2p\x\p\ + 2p 2 ^ 2 ]92) 
(since PjX 2 pj > 0,j — 1,2) 









> 


(- 








(j 




+4^0" sina 








> 


(- 








(j 



[•••]- 29fa| cos a| [(pi =F ^i) 2 - p\ - x\ + (p 2 =F x 2 ) 2 -p\- x\] 
[••■]+ 2Sa\ cos a I (p 2 + x\ + p 2 2 + x 2 2 ) + 43a sin a. (B.2) 



We shall prove below that the term inside square brackets in (B.2) satisfies the following 
estimate: 



2 1 i4 



[-e ta A + \(3\V{x)){-e- ia A + \(3\V(x)) + \a\ 2 \x 

|2 



> aM +pj) 2 + \0\ 2 V{x) 2 } + ^L|x| 4 - h. 



(B-3) 



for suitable constants ai,b\ > 0, independent of 7 G Vt and a in a compact subset of 

]0,7T[. 

Now, using (B.3) and setting A = (1 - ^ )%, 5 = \(l - $f ) and b 2 = (1 - ^ )&i - 
43 ; cr since, (B.2) can be bounded from below by: 

MpI + Pi? + \P\ 2 V(x) 2 ] + B\a\ 2 \x\ 4 - b 2 
+2^sa\ cos a\ (p 2 + x 2 + pi + x 2 ) 

Aa'Qo 2 + p 2 2 ) 2 + 29f<r| cos a| {p\ + p 2 2 ) - b 2 + 



b' 
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b' 

a'\B\\a\ 2 \x\ 4 + 29f<r| cosa|(:r 2 + x\) + - 



V(x) 2 + 

- a')(p? + pi? + B(l - a')\a\ 2 \x\ 4 - b' 



(B.4) 



Since the terms inside square brackets in (B.4) are positive for a suitable choice of the 
constants a', b' > 0, a' < 1, we finally obtain 

(-e ta A + ax 2 + \[3\V(x))(-e- ta A + ax 2 + \P\V(x)) 

> A(l - a')(pl + vlf + A\(3\ 2 V(x) 2 + B(l - a')\a\ 2 \x\ 4 - b' . (B.5) 

Now (B.l) follows from (B.5) with a = min(A(l — a'), B{1 — a')) and b — ^. In order 
to complete the proof of the lemma we need to prove (B.3). We have 

I 1 2 

(_ e -A + \(3\V(x))(-e- ia A + \(3\V(x)) + l -^-\x\ 4 

= (Pl+Pi) 2 + \P\ 2 V(x) 2 + \P\ cosa[(p 2 +p 2 2 )V(x) + V(x)(p\ + P \)] 

I 1 2 

+i\P\sma[{p 2 1 +p 2 2 ),V(x)]+ l -^-\x\ 4 

= {p\+pI? + \P\ 2 V(x) 2 ± oosalKp? +pl)V(x) + V{x){p\ + pj)} 

1 



-2 sin a 



x 2 {pixi + xipi) + xi(p 2 x 2 + X2P2) + ^{xjp2 +P2%1) 



= I cosa\(pi+p 2 ± \(3\V{x)Y + (1 - I cosa|)[(pi + p 2 ) 2 + \(3\ 2 V(x) 

X2{P\Xl + + Xi(p 2 X 2 + X2P2) + ^(^2 +P2X 2 2 ) 



-2 sin a I 



+ 



+ 



km 



.7; 



> (1 - I cosa|)[(p 2 + pi) 2 + \P\ 2 V(x) 2 ] + ^\x\ 4 + 2\/3\ sin a ( Pl - x^f 
+ (p 2 - x x x 2 ) 2 + ^(p 2 - x 2 2 ) 2 - [p\ + -p\ - 2x\x\ - ^x 4 

> (1 - I cosa\)[(pl + pi) 2 + \P\ 2 V(x) 2 } - 2\f3\ sin a (p\ + ^ + 2^2 + ^ 



KM 



\x\ 



(1- |cosa|)a 2 (p?+p^) 2 -2|^|sina^ + 7^2) + &i 
+ (1 - I cosa|)(l - a 2 )(p 2 + P2) 2 + (1 " I cosa|)|/5|V(x) 2 - h 



km 



|x| — sina(4a; 1 a; 2 + x 2 )- 



(B.6) 



Now, for a suitable choice of the constants < 02 < 1,63 > 0, the term in square 
brackets in (B.6) is positive and therefore (B.6) can be bounded from below by: 

(1 - I cosa|)(l - a 2 )(p 2 1 +p 2 2 ) 



KM 



+ (1- \cosa\)\(3\ 2 V(xY - b x + ^-\x\ 4 - \(3\sma(4xfx 2 2 + x 4 ). (B.7) 
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■|2 



Next notice that, for \(3\ sin a < we have 



12 

— \x\ A — A\(3\ sinax^xl — \(3\ sinax^ 



-x\+ ( _ |/5|sina] 4 + (kl 2 -4|/3|sina)x^ > 0. 



2 

Thus, we finally obtain 

(_ e -A + |/3|\/(x))(-e- itt A + |/3|F(a;)) + |a| 2 |x| 4 

> (1 - | cosa|)(l - a 2 )(pj +p 2 2 ) 2 + (1 - | cosa|)\/(:r) 2 

+ ^\A A ~h (B.8) 
which corresponds to (B.3) with a\ = (1 — | cosa|)(l — 02). 
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